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We explain the motivation and main results of our work in Ref. 1. Using the covariant 
formalism, we derive the equations of motion for adiabatic and entropy perturbations at 
third order in perturbation theory for cosmological models involving two scalar fields, 
and use these equations to calculate the trispectrum of ekpyrotic and cyclic models. The 
non-linearity parameters /jvl and g^L a-re found to combine to leave a very distinct 
observational imprint. 

Keywords: Early universe cosmology; Non-Gaussianity. 

1. Introduction 

The analysis of non-Gaussianity,^ i.e. the deviation from perfect Gaussian statis- 
tics, is becoming a strong discriminator of competing models of the early universe. 
Whilst most efforts until now have been devoted to the study of the bispectrum 
- the Fourier transform of the 3-point function of the curvature perturbation ^ - 
recent progress in observational cosmology is such that it may be possible in future 
to detect the trispectrum (4-point correlation function of C), if it is large enough. 
However, the calculational technology required to make the corresponding predic- 
tions is not yet in place. In particular, analytic formulae can be obtained from the 
5N formalism, the most widely used technique, only in very special cases and fur- 
thermore its numerical implementation is demanding. Here we report on a recent 
development: the derivation of the equations of motion to third order in perturbation 
theory, which are needed to calculate the non-Gaussianity in two-scalar-field cosmo- 
logical models. The resulting equations are generally applicable and well suited to 
numerical analyses, as demonstrated by our calculation of the trispectrum generated 
by ekpyrotic and cyclic cosmological models — scenarios in which the primordial 
fluctuations are created during a collapsing phase before the Big-Bang (see Ref. 3 
for a review). 



2. Multifield cosmological perturbations 

We consider cosmological models involving two fields with standard kinetic terms. 
In such models large non-Gaussianities can be generated because of the classical 
nonlinear evolution outside the horizon: in general multiple field models, the scalar 
perturbations can be decomposed into (instantaneous) adiabatic and entropy modes 
by projecting, respectively, parallel and perpendicular to the background trajectory 
in field space.* As the entropy fields are not submitted to the slow-roll requirements, 
they need not be almost free fields and can develop large nonlinearities. If they are 
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light enough to be quantum mechanically excited during inflation, they develop 
super-Hubble fluctuations that can be transferred, together with their nonlinear- 
ities, to the adiabatic mode or curvature perturbation, if there is a turn in fleld 
space (or even more generically with non standard kinetic terms^). The conversion 
from entropy to curvature perturbations can also itself be nonlinear. This results 
in possibly large non-Gaussianities of the so-called local form, characterized in the 
simplest models by two non-linearity parameters /j^£ and qnl via an expansion of 
the curvature perturbation Q in terms of its linear, gaussian part Cz,, 

c = a + If^ia + ^9NLa- (1) 

The determination of (resp. gNh) requires the derivation of the coupled equa- 
tions of motion for adiabatic and entropy perturbations on large scales to second 
order (resp. third order) in perturbation theory. This rather intricate task can be 
fulfllled by means of the covariant formalism, in which one easily derives exact fully 
nonlinear equations for covcictors that can then be expanded into perturbation the- 
ory. Langlois and Vernizzi used this method to derive the relevant equations to 
second order in Ref. 6 (see also Ref. 7), and their results were subsequently ap- 
plied to the determination of the bispectrum (/^£) in ekpyrotic and cyclic models 
in Refs. 8,9. In Ref. 1, Jean-Luc Lehners and I have extended these two works to 
third order in perturbation theory. It is worth emphasizing that our analysis is quite 
general: we did not specify an underlying potential and do not use approximations 
other than the large scale limit. In particular, although we have applied the formal- 
ism to a particular scenario of a contracting universe, it is equally well suited to the 
study of cosmological perturbations generated during a period of inflation. 

3. The ekpyrotic trispectrum 

In ekpyrotic scenarios, the curvature perturbation is generated through the entropic 

mechanism explained above. Broadly speaking, there are two limiting cases that are 
of special interest: the first is where the bending of the trajectory occurs after the 
ekpyrotic phase, during the approach to the big crunch. In this case, the kinetic 
energy of the scalar fields is the dominant contribution to the total energy density 
while the conversion takes place: we call this case "kinetic conversion" . The other 
limiting case is where the conversion occurs during the ekpyrotic phase — "ekpyrotic 
conversion" . During the ekpyrotic phase, we adopt the following parametrization of 
the potential: 

= -Foe^'^[l + es' + |e^/^s3 + ^eh' + •••], (2) 

where a denotes the adiabatic direction, s denotes the "entropy" direction, we 
expect K3,K4 ^ 0{1) and where e ~ O(IO^) is related to the ekpyrotic equation of 
state Wek via e = 3(1 -|- Wek)/^- 

Then, for kinetic conversions lasting of the order of one e-fold of contraction of 
the scale factor, we find the following approximate fitting formula (for which there 
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is an analytical understanding^'^) for the third order non-linearity parameter qml 
(and we include for completeness®^'') : 



Note that when ^3 and Ki are small, qnl is always negative. Hence, even though 
/'y£ is small in that case, g^L is negative and typically of order a few thousand, 
so that any accidental degeneracy at the level of between simple inflationary 
and cyclic models is very likely to be broken at the level of the trispectrum. More 
generally, unless |/^£| turns out to be quite large, one would typically expect gjvi, 
to be negative, as obtaining a positive qnl in that case would require unnaturally 
large values of K4. 

For ekpyrotic conversions, the SN formalism is easily applicable and gives ap- 
proximate formulae for the non-linearity parameters that we were able to confirm to 
a good accuracy by numerically solving the equations of motion up to third order. 
In that case /^£ is always large and negative while gNL is always positive. The 
current limits on fj^l^^ are therefore already putting some strain on this particular 
mode of conversion. 

4. Conclusion 

Motivated by huge observational perspectives in coming years concerning non- 
Gaussianity measurements, we have developed in Ref. 1 the necessary tools to 
compute the trispectrum generated in two-field cosmological models. Our results 
are generally applicable and well suited for numerical analyses and we used them to 
calculate the trispectrum generated in ekpyrotic and cyclic models. The combined 
consideration of the bi- and trispectrum in that case results in a distinct observa- 
tional imprint that should enable one to select or rule these models on observational 
grounds in the foreseeable future. 
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